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Abstract 

We analyze the power output of a quantum dot machine coupled to two electronic reservoirs via thermoelectric contacts, and to 
two thermal reservoirs - one hot and one cold. This machine is a nanoscale analogue of a conventional thermocouple heat-engine, 
in which the active region being heated is unavoidably also exchanging heat with its cold environment. Heat exchange between 
the dot and the thermal reservoirs is treated as a capacitive coupling to electronic fluctuations in localized levels, modeled as two 
additional quantum dots. The resulting multiple-dot setup is described using a master equation approach. We observe an “exotic” 
power generation, which remains finite even when the heat absorbed from the thermal reservoirs is zero (in other words the heat 
coming from the hot reservoir all escapes into the cold environment). This effect can be understood in terms of a non-local effect 
in which the heat flow from heat source to the cold environment generates power via a mechanism which we refer to as Coulomb 
heat drag. It relies on the fact that there is no relaxation in the quantum dot system, so electrons within it have a non-thermal 
energy distribution. More poetically, one can say that we find a spatial separation of the first-law of thermodynamics (heat to work 
conversion) from the second-law of thermodynamics (generation of entropy). We present circumstances in which this non-thermal 
system can generate more power than any conventional macroscopic thermocouple (with local thermalization), even when the latter 
works with Carnot efficiency. 


Highlights 


• Model of a quantum-dot energy-harvester; a quantum analogue of a thermocouple. 

• Two thermal reservoirs (heat source and cold environment) and two electric contacts. 

• Non-local Coulomb heat drag effects, non-locality of laws of thermodynamics. 

• Finite power output without absorbing heat from the thermal reservoirs. 

• Can beat a Camot-efficient conventional thermocouple under equivalent conditions. 
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Dedicated to Markus Biittiker: In addition to his human 
qualities, we remember Markus fondly for the inspiring discus¬ 
sions we had. We think that this work would have led to another 
lively and enjoyable debate. 

1. Introduction 

There is great current interest in quantum and nanoscale sys¬ 
tems that convert heat into electric current mB- The simplest 
such systems are those that exhibit a thermoelectric effect EUS: 
they do the conversion in a steady-state (DC) manner, and so 
avoid the need for pumping cycles relying on time-dependent 
couplings. Quantum dots are particularly promising in this re¬ 
spect, and various applications have been proposed and (at least 
partially) realized experimentally, including thermoelectric en¬ 
gines llbi EOl . refrigerators ETT - ESl . thermal rectifiers Il26ll27l . 
and hybrid refrigerator power-sources ||28]. In addition, the 
simplicity of these thermoelectric quantum dots makes them 
ideal model systems for the study of quantum thermodynam¬ 
ics. 
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The most developed theory of non-equilibrium thermody¬ 
namics, known as irreversible thermodynamics Il29l . assumes 
that the shortest lengthscale in the system is that on which the 
particles thermalize (the inelastic scattering length). Then, any 
system coupled between two reservoirs at very different temper¬ 
atures has a well-defined local temperature everywhere within 
it. However, as we reduce the size of the circuit elements into 
the nanoscale regime, this ceases to be the case. A system that is 
much smaller than the electron thermalization length will be in 
a non-thermal state whenever it is coupled to reservoirs that are 
at significantly different temperatures or electrochemical poten¬ 
tials. In other words, the distribution of the excitations in the 
system cannot be described in terms of a thermal distribution. 

In general, we wish to answer the question; what new physics 
can emerge in a quantum system operating in this non-thermal 
regime? This work presents a first response to this question, in 
the context of the thermoelectric quantum dot system sketched 
in Fig. [^a). This can be considered as a miniaturized ver¬ 
sion of the usual macroscopic thermocouple power generator, 
in which the macroscopic metal reservoir between the two ther- 
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(a) Quantum dot thermocouple 



(b) The situation shown in (a) drawn as a circuit 
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Figure 1: (a) A sketch of a quantum dot thermocouple in which the electrons 
in the central part of the thermocouple (the quantum dot) exchange heat with its 
cold environment (Reservoir C) as well as the heat source (Reservoir H). The 
couplings of the dot to reservoirs L and R have an energy dependence chosen to 
ensure that each one has the opposite thermoelectric response. As a result they 
form a thermocouple, which can be used to generate electrical power. The load 
is taken to be a device that converts the electrical power into some other form of 
useful work (it could be a motor converting the electrical work into mechanical 
work), (b) A circuit picture of the set-up in (a), in which the exchange of energy 
with the heat source and the environment are modeled by a capacitive coupling 
to electronic fluctuations in those reservoirs, which we take to be thermally 
activated hopping of electrons between the bulk and a localized state, indicated 
by the two upper circles (H and C). Thus each of the three circles represents a 
two-state system, with occupations 0 or 1. 


moelectrics has been replaced by the quantum dot M. Our sys¬ 
tem is similar to three-terminal energy harvesters considered 
in Refs. lH) 1^ . which separate the conductor from a heats 
source with which it exchanges energy but no particles. See also 
Refs. ElllIHSa, which consider related models with bosonic 
heat sources. However, instead of the dot being coupled only 
to the heat source (reservoir H), we consider the case when it 
is also coupled to the cold environment (reservoir C). Since the 
environment coupling is never negligible, we argue that this is 
the generic case. It is certainly the case in conventional thermo¬ 
couples, where it is well known that the region between the two 


thermoelectrics is not as hot as the heat-source, because it also 
exchanges heat with the cold environment. With this in mind, 
the simplest case would be when the temperature of the cold en¬ 
vironment C is the same as the temperature To of the electrical 
circuit in which the thermocouple is inserted (reservoirs L and 
R). However, we consider the more general case where Tq is 
different from Tq, for example due to Joule heating in the wires 
or load {Tq> Tq). Following what is known about conventional 
thermocouples, one would expect the energy available to pro¬ 
duce electricity to be equal to the heat that flows in from reser¬ 
voir H minus the heat that is lost as it flows out into reservoir 
C. We call this quantity 7in, and study how the power generated 
by a quantum dot thermocouple, Rgen, depends on it. 

In this work, we show that the quantum-dot thermocouple 
can generate power even when the total heat absorbed from the 
thermal reservoirs (H and C) is zero (7in = 0). Instead, the dot 
extracts heat from the electronic reservoirs (L and R) and con¬ 
verts it into electrical power. It might appear that this “exotic” 
power generation violates the laws of thermodynamics (since 
the thermal reservoirs provide no energy to allow the dot to con¬ 
vert heat into work), but we show that this is not the case. It can 
be explained in terms of a non-local heat drag effect, in which 
the heat flow from H to C can induce heat and charge currents 
in the circuit (i.e. through the thermoelectric contacts to reser¬ 
voirs R and L) even when 7in = 0. We show that this exotic 
power generation only arises because no thermalization occurs 
within the quantum dot, meaning that the dot is maintained in a 
non-thermal state by its contact with the hot and cold reservoirs. 
As such it has no analogue in conventional thermocouples. We 
argue that as a result a quantum-dot thermocouple can achieve 
power outputs larger than any conventional thermocouple (even 
one working with Carnot efficiency) over a broad range of pa¬ 
rameters. This is not because our device has an efficiency higher 
than Carnot efficiency (this is forbidden by the laws of thermo¬ 
dynamics), but because it can also extract useful work from the 
non-local drag effect. This is something a conventional thermo¬ 
couple cannot do, irrespective of how efficient it is. 

This work is organized as follows. In Sec. we introduce 
our model of the quantum dot thermocouple. Sec. [^discusses 
the basic mechanisms behind the rectification of heat fluctua¬ 
tions in our device, in the simpler case where it generates no 
power. Sec. gives results for the power generation, which 
Sec. 1^ explains in terms of non-local heat drag, shows that it 
obeys the laws of thermodynamics, and discusses potential ex¬ 
perimental implementation. Sec.j^shows that the effect is sup¬ 
pressed if relaxation processes cause the state of dot M to be¬ 
come thermal. Finally, in Sec. [7] we compare the power output 
of the quantum dot thermocouple with that of a conventional 
thermocouple for a broad range of parameters. 

2. Model of the quantum-dot thermocouple 

The quantum dot M, shown in Fig.[TJa), can exchange energy 
(but not charge) with a hot (H) and a cold (C) thermal reservoir. 
In the following we are interested in the resulting amount of 
electrical power that this quantum dot can produce by driving 
an electric current between the left (L) and right (R) electron 
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reservoirs, to which it is tunnel-coupled. Importantly, while the 
reservoirs are characterized by temperatures, the quantum dot 
is generally in a non-thennal state. This is a crucial point of our 
paper and an important difference to conventional macroscopic 
thermocouples, which are large enough such that every element 
of the thermocouple is able to thermalize. 

The energy exchange of the quantum dot M with the thermal 
reservoirs H and C occurs via a capacitative coupling to thermal 
fluctuations in these reservoirs. They are modeled as thermally 
activated hopping of electrons between the bulk of each of the 
two thermal reservoirs and a localized state on their surface. 
The localized states can be thought of as quantum dots, which 
enables us to recast the problem as the triple-dot circuit shown 
in Fig.[2b). Here, the localized states associated with reservoirs 
H and C are the upper red and blue dots. Apart from thermo¬ 
electric purposes, similar models have been useful for relax¬ 
ation times detection |(38], Maxwell demon physics ll^ l40l , 
and the operation of stochastic logic gates HTI . 

We therefore now consider the three capacitively coupled 
quantum dots M, H, and C, in contact with four electronic reser¬ 
voirs H, C, L, and R, as shown in Fig.[TJb). We assume that the 
on-site Coulomb repulsion is very large, such that doubly occu¬ 
pation of any dot can be neglected. For simplicity, we neglect 
the spin degree of freedom of the electrons, which, in this case, 
does not alter the working principle of the device, see also 0 
(note however, that in general the spin might play an impor¬ 
tant role in quantum-dot based heat engines ||42l). The isolated 
system of three capacitively coupled quantum dots can then be 
described by the model Hamiltonian 

■K = ^ eihi+ ^ Uijhihj, (1) 

i=M,H,C ;,j=M,H,C 

with the number operator n, counting the number of electrons 
(0 or 1) on each of the three dots, i - M,H,C. Here, e, are the 
energies for single occupation of the different dots and Uij rep¬ 
resents the Coulomb charging energy that needs to be paid to 
have both dots i and j occupied simultaneously. They play an 
essential role, as they mediate the energy exchange between the 
different dots. The eigenstates of this Hamiltonian are given by 
the eight states (m, h, c) with eigenenergies where the la¬ 
bels m, h, c can take the values 0 (empty) and 1 (full) depending 
on the occupation of the respective dots M, H, and C. 

Transitions between eigenstates of the triple-dot system can 
occur due to the tunnel coupling of these dots to the electronic 
reservoirs. These reservoirs are kept at temperatures £„ and 
electrochemical potentials fia- Here we take the temperature 
of the hot reservoir to be always larger than the temperature 
of the cold reservoir, Th > Tq- The left and right reservoirs 
are assumed to be at the same temperature, Tl = Tr = Tq 
for Tq < Th. In this paper, we are interested in the power 
that the quantum-dot thermocouple can generate, by driving an 
electrical current against a potential difference (due to a load). 
A simple way of implementing this setting is the one where 
we take the electrochemical potential of the right lead to be 
fiR - eV > 0 (where e < 0 is the electron charge) with respect 
to the equal electrochemical potentials of the other reservoirs, 
which we take as the reference energy here, jUh = yuc = A<l = 0. 


The coupling between the quantum dots and the reservoirs, 
as depicted in Fig. [^b), is characterized by the tunnel cou¬ 
pling strengths, ja- We focus on a situation, where the dots 
are weakly coupled to the reservoirs, ja ^ k^Ta- Importantly, 
since we are interested in the quantum dot M behaving as a 
thermoelectric, we require that -/l and yR depend on energy 
Because of the capacitive coupling between the dots, this means 
in turn that they depend on the occupation h and c of the dots 
H and C, yielding yL = 7 Lih,c) and yR = yR(/!,c). A strong 
energy-dependence of the tunnel-coupling can be achieved for 
example by coupling the quantum dot to the reservoirs via other 
resonant levels EBS- This is not required for the tunnel¬ 
coupling between the thermal baths and the nearby quantum 
dots. We therefore assume yn and yc to be energy-independent. 

We are interested in the steady-state response of the system to 
different temperature in the heat baths, in the presence of a pos¬ 
sible voltage drop across the reservoirs L and R, which results 
in power being generated across the load. In order to calculate 
the DC charge and heat currents, the occupation probabilities 
Pmhc of the states (m, h, c) of the triple-dot system are required. 
In the weak-coupling regime, these probabilities are given by 
the solution of the steady-state master equation llT5ll44ll . 

0 = X , (2) 

m' ,h' ,c' 


with Yjmhc Pmhc - 1 ■ All non-diagonal elements of the transition 
matrix are determined by 


m; 


mhc 




h'c 
^-1 mhc 


Q'=H,C,L,R 


( 3 ) 


namely by the sum over the rates of transition of the three-dot 
system from state m'h'c' to state mhc due to electron tunnelling 
into or out of reservoir a, [Fq,] . The latter are given by 



), (4a) 

r -1 m'h'c' / nil 

mhc ~ ^m,\-m'5h,h'^c,c' 7R{h, c)fR {j^mhc 

-(m-m')eV^, 


(4b) 

~ ^>n,m'^h,l-h'Sc,c' 7iifR[^mhc )’ 

(4c) 

M2"- <m'5hM6n^c' yc/c (Alf'), 

(4d) 


where is a Kronecker delta, 

and fa{x) - 1/(1 + pej^j function correspond¬ 

ing to lead a. The diagonal elements in Eq. Q are directly 
found by using the condition that every column of the tran¬ 
sition matrix sums to zero to fulfill probability conservation, 

jiyfmhc _ _ v 

For any heat-engine, the relevant quantities are the generated 
power (a charge current flowing against a potential difference) 
and the heat flow into the device. We therefore consider the 
charge and heat currents into the dots from reservoir a, 


la - el^ (5a) 

- Pala- (5b) 
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Here the charge current is trivially related to the particle cur¬ 
rent, I^, by multiplication with the electronic charge e, while 
the heat current is given by the rate with which energy changes, 
/f, with respect to the energy that particles at the electrochem¬ 
ical potential would carry in the respective reservoir. Using the 
Clausius dehnition of entropy, we see that the rate of change of 
entropy in reservoir a is connected to the heat current by 


(d5 Jdf) = -JJT,. (6) 

The particle currents and energy currents [^can be written using 

[ -1 m'h'c' 

Tq. and the occupation probabilities 

-I mhc 

Pmhc of the state (m, h, c) with the respective energy Emhc and 
occupation number N,nhc - m + h + c, 

j j i^mhc Pm'h'c' ? (7) 

mhc m'h'c' 

^ ^ Pm'h'c' ■ (8) 

mhc m'h'c' 


The particle and energy currents in reservoir a are dehned to be 
positive when they are directed flowing into the triple-dot sys¬ 
tem. Consequently, the charge current in reservoir a, Eq. (5ai, 
is positive when electrons flow into the reservoir, while the heat 
current in reservoir a, Eq. (5b i, is positive when flowing into the 
triple-dot system. Particle and energy currents are conserved, so 
summing a over all reservoirs gives 


0 = (9) 

a a 


Charge current is conserved (so /r = -/l)- Heat current is 
not conserved since which corresponds to 

the hrst law of thermodynamics. The power generated by the 
quantum-dot thermocouple is 


Eor strong t/nc^ the states (0,1,1) and (1,1,1) drop out of the 
problem. In this limit, the charge fluctuations that lead to trans¬ 
port through the conductor can be separated in two cycles. The 
first cycle is 

(0,0,0) ^ (1,0,0) ^ (1,1,0) ^ (0,1,0) ^ (0,0,0), (11a) 

and involves an energy exchange I/mh between the system and 
the hot reservoirs. The second cycle is 

(0,0,0) ^ (1,0,0) ^ (1,0,1) ^ (0,0,1) ^ (0,0,0), (11b) 


and involves an energy exchange Uuc between the system and 
the cold reservoir. When the tunneling events in these se¬ 
quences involve an electron that tunnels from reservoir L into 
the dot, and thereafter into reservoir R (or viceversa), it con¬ 
tributes to transport. The sum of all tunneling processes results 
in the generation of a finite electric current by the mere con¬ 
version of heat, only if a preferred direction is defined both by 
broken electron-hole and left-right symmetries 0. We impose 
these conditions on the corresponding tunneling rates, requiring 
that the asymmetries, 

yL(l,0)yR(0,0)-yL(0,0)yR(l,0) 

“ [yL(0,0) + yR(0,0)][yL(l,0) + yR(l,0)] ’ 

and Ac (obtained by replacing yQ.(l, 0) ^ yff(0,1) in the previ¬ 
ous expression) are finite. 

Solving the master equation Q, we get analytical expres¬ 
sions for the charge and heat currents. Despite these expres¬ 
sions being long and cumbersome (even for V - 0), we And a 
remarkably simple relation between them: 


/R(y=0) = e 


Ah 


7h + 



(13) 


It is instructive at this point to introduce and 7trans, where 


Egen = -PrI^ = -VIr . (10) 


■^in - + -^C . 


(14a) 


These quantities are evaluated using Eqs. Q, and we will dis¬ 
cuss them in detail in the following sections. Eor now, we do 
not consider the probe reservoir P, shown in Pig.[^ however we 
will discuss it in Sections |2and|7] 

3. Results in the absence of a bias voltage 

We start by analysing the behaviour of the quantum dot ther¬ 
mocouple in Pig. [funder the following two assumptions, (i) 
The load’s resistance is zero, so the bias V remains zero when a 
current is flowing, (ii) The capacitive coupling between dot H 
and C is strong enough that Uhc is very much bigger than all 
temperatures, so the probability that these two dots are occu¬ 
pied at the same time becomes vanishingly small. We will relax 
both these assumption in Section]^ 


* These currents are time-averaged, and do not include the noise that comes 
from short-time fluctuations. This noise has a negligible elfect on the long¬ 
time (steady-state) conversion of heat into work, see for example Section III of 

Ref. (TT). 


is the total heat flowing into dot M from the thermal reservoirs 
(reservoirs H and C), and 

Arans ~ “ ■fc ■ (l^h) 


is the heat that is transferred from reservoir H to C. Note that a 
proportion of 7trans may transit through dot M on its way from 
reservoir H to C, even when none of it enters the electrical cir¬ 
cuit (i.e. reservoirs L, R). Rewriting Eq. ( [T3] l in terms of 7in and 
■/trans, One sees that /R(y=0) contains one part that depends on 
Jin as one might expect, but it also contains a term that depends 
on 7trans, even though that energy does not go into the circuit. 

If the temperature and coupling of the hot and cold reservoirs 
are such that Jin = 0, one still generates a charge current [which 
we here denote by I* = I{Jin - 0)] in the circuit equal to 


r(y=o) 


g / Ah 
2 \ Umh 



(15) 


Whenever the charging energies Umh and Umc are not equal, 
there is no symmetry reason for I* to be zero, so a typical sys¬ 
tem will exhibit a charge current in the circuit even though there 
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Figure 2: Upper plot: Power generated by the quantum-dot thermocouple 
(solid curve) as a function of Tq compared to the maximum power T/camot^fin 
(dashed line) any conventional thermocouple could generate (see Sec. [TJ. 
Lower plot: Energy currents into the system as a function of TqITq. For both 
panels, the system parameters are given by the temperatures k^Tuly = 16, 
k^Toly = 8, the dot energy levels 6^ = = 0^ the charging energies 

t^MH/y = 3, Uucly = 6, Uwdy = lOO, and the coupling constants yn/y = 1. 
ycjy = 0.5, while yL and yR are given in Eq. (El¬ 


is no energy absorbed from the thermal baths (7in = 0). In con¬ 
trast, in the special case of I/mh = Uuc, one can show that 
Ah = Ac, which means that /*(y=0) = 0. To show this, one 
notes that if Umh = Umc, then the tunneling rates into dot M 
are not sensitive to which of the upper dot is occupied (H or 
C), so 7 q,( 1,0) = yQ-(0,1), which implies Ah = Ac- It is thus 
essential that the two baths act differently on the system. We 
achieve this by taking different capacitive couplings. 


4. Results for finite power generation 


For a finite power output, the machine must drive a current 
against a finite potential difference V. This occurs when the 
load’s resistance is non-zero. We also now consider the general 
case in which t/nc is finite, so the dot dynamics involves more 
cycles than just those in Eqs. (Ill, which can explore all eight 
possible states of the three-dot system. Since in this general 
situation at any V, the algebraic solution of the master equa¬ 
tion, Eq. and of the resulting charge and heat currents are 
extremely cumbersome, we refrain from presenting them here. 
Instead, in Eigs. m we plot the power output generated by our 
quantum-dot thermocouple as a function of the temperature of 
the cold reservoir, Tq, for various specific system parameters. 
Eor comparison, we also show the energy currents flowing from 
reservoir L and R into the dot M, and /^, as well as the heat 
current flowing into the system from the hot and cold reservoirs, 
An- 


These plots show the possibility of generating finite power 
when the energy flow into the system from the hot and the 
cold reservoir vanishes. Importantly, the effect is shown to 
be generic, since the variation of different system parameters 


does not alter the basic features. Eor the plots in this paper, we 
give all rates and energies in terms of a rate y. As a result the 
power generated Pgen and all energy currrents Jo are in units of 
y^ (we take Planck’s contant H = 1). In all cases, we fix the 
temperatures of the hot reservoir and the load, k^Tuly = 16, 
ksToly - 8, and vary k^Tc. In different plots we take differ¬ 
ent tunneling rates, {yo}, but always take them to be smaller 
than the temperatures, which ensures the applicability of the 
master equation approach introduced in the previous section. 
In each plot, for each value of Tq, we find the bias voltage 
V that maximizes the power generated. We call this optimal 
bias Vopt(Tc)- We then plot the power generated and the energy 
flows at y = yopt(7’c) for each Tq. 

The plot in Eigs.andare for 

'Yi.(h,c)ly = Sh,o6c,o, (16a) 

7 Rih,c)ly = l-dh,o^c,o- (16b) 


This choice of yi^{h, c) and ys.(h, c), ensures that when dot H 
and C are empty dot M is decoupled from reservoir R, but if ei¬ 
ther dot H or C is occupied then dot M is decoupled from reser¬ 
voir L. This means that the asymmetry parameters in Eq. © 
are maximal, with Ah = Ac = -1. All the other parameters are 
given in the figure captions. As mentioned before, our model 
contains too many parameters to present a systematic study of 
the full parameter space, however the features we find appear 
to be generic. The following two paragraphs highlight this by 
discussing central points about Eigs.|^and[^ 

In Eig.|^we take a very large Uhc (much larger than ^bTh), 
which means that the dots H and C are so strongly capaci- 
tively coupled that there is an extremely small probability of 
them being both occupied at the same time (similar to the limit 
f/nc ^ “ taken in Sec. |^. The power generated, Pgen, van¬ 
ishes quadratically at a certain value of Tc. This parabolic be¬ 
haviour of Pgen = -JrV, occurs because the current /r changes 
sign (going linearly through zero), and thus the bias generated 
also changes sign at this point (y oc /r for small /r). However, 
we clearly see that Pgen vanishes at smaller Tq than Ji„; so Pgen 
is finite when Ji^ - 0 (marked by the filled circle in Pig. |^. 
The power generated at this point is purely an effect of the heat 
flow 7trans^ which passes from reservoir H to C without being 
absorbed by the thermocouple circuit. We call this the exotic 
power generation. When 7in = 0, one has - 0, thus 

finite Pgen implies that there is a net heat current out of reser¬ 
voirs L and R, meaning that the total entropy of those reservoirs 
is dropping. Sec. 5.2 explains why this is not a violation of 


the second law of thermodynamics|^ We do not show 7tran.s in 
Pig. HI as it is significantly larger than the other heat flows, it is 
a positive and monotonically decreasing function of Pc/Po (it 
goes from 0.14 at Pc/Po = 0.2 to 0.03 at Pc/Po = 1-2). 

The plots in Pig.|^show the same basic features as in Pig.|^ 
despite the following differences in the choice of parameters. 


^Fig.j^has = 0. however this has no physical consequences as is 
gauge-dependent (it depends on our choice of E = 0). In contrast, the other 
quantities considered above, such as (//^+/^), P and are gauge-independent 
and so have physical meaning. 
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Figure 3: Upper plot: Power generated by the quantum-dot thermocouple 
(solid curve) as a function of Tq compared to the maximum power T/camot'^in 
(dashed line) any conventional thermocouple could generate (see Sec. [TJ. 
Lower plot: Energy currents into the system as a function of TqITq. For both 
panels, the system parameters are given by the temperatures = 16, 

^bT’o/'X = 8, the dot energy levels 6 m = = 7y, the charging energies 

Umh/t = 5, Vudy = 6, Uwcly = 3, and the coupling constants ynly = 0-5, 
jQjy = 1, while 7 l and 7 r are given in Eq. i|16k. The bias V = Vopt(Tc) is the 
one that maximizes the power generated for each TqITq. 

We no longer take [/hc ^ be effectively infinite (i.e. much larger 
than kBTii). We take the energy-levels of the dot to be above the 
chemical potentials of the reservoirs. The coupling to reservoirs 
L and R are still given by Eq. ensuring the optimal asym¬ 
metry for the thermoelectric effect. The decreased capacitive 
coupling Uiic (of the order of temperature) allows all the states 
of the three dots to play an active role in the dynamics (unlike in 
the limiting case Uuc °° considered in Sec.[^. This however 
appears to have little qualitative effect on the physics (one can 
even take Uhc - 0 without considerably changing the features 
that interest us). However, the shift of the energy levels with 
respect to the electrochemical potentials is directly responsible 
for the only qualitative difference between the results in Figs.l^ 
and 1^ namely that Tin now vanishes at two points in Fig. 
This can be explained as follows. When Tc is large, then Tin is 
positive, but as one reduces Tc, the heat flow Tin goes to zero 
and becomes negative. However, as one reduces Tq even further 
k^Tc becomes significantly less than ec, and so the dynamics of 
dot C start to freeze out (thermal fluctuations of dot C become 
extremely rare, as ec “ A'c ^ k^Tc). Thus, when Tq is suffi¬ 
ciently small, reservoir C becomes effectively decoupled from 
the rest of the system, the only heat flow is that from reservoir 
H to reservoirs L and R, and becomes positive again. The 
power generated remains positive at all Tq because of the exotic 
contribution coming from Ttrans- Another consequence of the 
off-resonant level positions chosen here, is that both the energy 
current out of the left and the right reservoir are in general non¬ 
zero. As required from energy conservation, they are however 
opposite, = -/^, at the point of vanishing 7in. Again 7trans 
is positive and significantly larger than the other heat flows for 



Figure 4: Upper plot: Power generated by the quantum-dot thermocouple 
(solid curve) as a function of Tq compared to the maximum power //Camot'fin 
(dashed line) any conventional thermocouple could generate (see Sec. |7|. 
Lower plot: Energy currents into the system as a function of TcITq. For both 
panels, the system parameters are given by the temperatures k^Tnly = 16, 
k-hToly = 8, the dot energy levels 6 m = = 6y, the charging ener¬ 

gies Umh/t = k Umc/t = 12, and U^dy = 100. and the coupling constants 
yB = yc = y. while yL and yR are given in Eq. jl7|. The bias V = Vopt(Tc) is 
the one that maximizes the power generated for each TqITq. 

all TcITq, as we show in the inset of Fig. It drops for small 
Tc, because of the above mentioned freezing of the dynamics 
of dot C. 

The plots in Fig. are for an imperfect thermocouple, where 
dot M is coupled to both reservoir L and reservoir R for all 
occupations of the of upper dots {h and c). We replace Eq. 
with 

yc(.h,c)ly = (2d;,,o<Jc,o)/3, (17a) 

yRih,c)ly = 0-6f,fiScfi)l3. (17b) 

This means that the coupling strength is only varying by 1/3 as 
h and c change and the asymmetry parameters in Eq. are 
Ac = Ah = -0.2. The other parameters chosen for this plot 
are given in the figure caption. The figure shows that the exotic 
power generation still exists for this imperfect pump, although 
its magnitude is greatly reduced (just as the traditional power 
generation is greatly reduced). 

5. Interpretation 

5.7. Coulomb heat drag 

The charge current we found in Sec. for 7in = 0 is due 
to a form of Coulomb heat drag, in which the heat flow from 
reservoir H to reservoir C, drags a heat and charge currents from 
reservoir L to R. This is possible due to the capacitive coupling 
between the quantum dots. However, in other drag effects |l8l 
|45]| . there is a transfer of energy associated with one current 
dragging the other, which corresponds to finite 7in. Hence, this 
drag effect without energy transfer is extremely unusual. 
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Figure 5: The entropy generated in each reservoir (main plot) and /tram (inset) 
for the system in Fig. The entropy generated in reservoirs L, R and FI are 
small in magnitude and their sum is negative, however the entropy generated in 
reservoir C is large and positive. Thus, the second law of thermodynamics is 
not violated. We see that the rate of entropy production broadly follows /u-ans 
(inset). We do not show plots of entropy production and /trans for the system in 
Fig.|^ as the main features are generic to all of them. 


This heat drag effect can occur for any Tc < Th irrespective 
of the value of TqITq, since it only requires a heat flow from 
reservoir H to reservoir C. The quantum dot thermocouple then 
produces more work than a classical one (see the discussion in 
Sec.7). Remarkably, this does not only occur in the configu¬ 
rations where Tin a; 0, but also in a relatively wide parameter 
range around this points, including the case Tc <: Tq. 

5.2. Non-locality and entropy production 

Despite the unusualness of the Coulomb drag without energy 
transfer discussed above, the situation is not paradoxical as long 
as y = 0. The case V = 0 corresponds to a load with zero re¬ 
sistance, so nothing opposes the charge flow, and so no work is 
required to generate a finite charge current. In contrast, as soon 
as the load’s resistance is finite, the charge current I must flow 
against a finite voltage bias, V, then work must be performed 
on the load. It would be natural to assume that this work must 
be provided by the thermal reservoirs, i.e. that the charge cur¬ 
rent I can only flow if Tin 0. However, something completely 
different can be observed in the results in Sec.|4| there we show 
that finite power can be generated even when Ji„ - 0. We call 
this exotic power generation. We observe that the power gener¬ 
ated is equal to the heat absorbed by dot M from reservoirs L 
and R, as there is no heat flow into dot M from anywhere else 
when Tin = 0. Thus the dot is sucking heat from these reser¬ 
voirs (which are both at the same temperature), to turn it into 
useful work. This obeys the first law of thermodynamics (en¬ 
ergy conservation), but - if we ignore the presence of reservoirs 
H and C - it appears to violate the second law, since it reduces 
the entropy of reservoirs L and R. 

However, if we now turn to reservoirs H and C, we see that 
there is a heat flow from hot to cold, and this increases the 
total entropy involved. An elegant way to verify this is the 
prove in Ref. ll46l that any such master equation respects the 
laws of thermodynamics, thus we know that the four reser¬ 
voirs and three dots together do produce entropy. What is 
more, one can rewrite the total rate of entropy production. 


(d5pr/df) = 2;„(d5„/df), as 


d5pr 

J in 1 

- 

1 

1 \ , '^trans , 

f 1 

_ M 

P 

-'gen 

dr 

2 ' 

[to 

Th 

Tc/' 2 

ITc 

Th/ 

To 


Now, since we know that (d5pr/df) > 0 for our system, see also 
Fig. 0 one knows that 


P < :^ 
^ gen ^ 2 


n 

Th 


To \ -/trans / Tq 

T^)^^\T~c 



(19) 


Hence, the second law allows the power generated, Pgen, to 
be finite and positive, even if Tin = 0, as long as Ttrans is fi¬ 
nite and positive. The entropy increase in reservoir C compen¬ 
sates for the reduction of entropy in the other reservoirs (see 
Fig.ig. Thus the apparent paradox is resolved by the fact the 
power generation has a non-local nature. The entropy genera¬ 
tion caused by the heat flow between one pair of reservoirs (H 
and C) enables the conversion of heat into work at another pair 
of reservoirs (L and R), even though there is no heat (or energy) 
flow between these two pairs of reservoirs and even though they 
can be arbitrarily far apart. 

A “poetic” way of expressing this non-locality is to say that 
there is a spatial separation of first and second laws of thermo¬ 
dynamics. The thermocouple manifests the first law by con¬ 
verting heat in reservoirs L and R into work. At the same time, 
some distance away, entropy is being generated by the flow 
of heat from reservoirs H and C, manifesting the second law. 
While we have shown that this behaviour is not paradoxical, it 
certainly is exotic! 


5.3. Towards experimental observation 

Transport through capacitively coupled quantum dots in mul¬ 
titerminal conductors has been measured IgTHlOl. The effect of 
large charging energies on the thermoelectric response of small 
tunnel junctions l34ll and quantum dots ll35]l is well-known ex¬ 
perimentally OISTJ. Our proposed device is close to those 
very recently realized in experiments lfT8l420l . in which a quan¬ 
tum dot generates a charge current between two reservoirs (L 
and R) by the rectification of thermal fluctuations in a third 
capacitatively-coupled reservoir (the heat source). In those sys¬ 
tems, it is likely that the dot with tunnel-coupling to reservoirs 
L and R also has at least a weak capacitative coupling to thermal 
fluctuations in its cold environment (including reservoirs L and 
R). Thus, these experimental systems have a strong similarity 
with the four-terminal set-up that we discuss. 

However, to unambiguously show the exotic power genera¬ 
tion experimentally, it would be desirable to have a more con¬ 
trolled coupling to the cold environment. This could be done by 
capacitatively coupling a third dot to the systems in Refs. IfTSl - 
|20l, which we believe is achievable with the technologies used 
in those works. More challenging will be to show experimen¬ 
tally that Tin = 0, for this it will be crucial to measure the heat 
currents accurately. For this, the capacitive coupling to the heat 
source and sink permits the different charge fluctuations to be 
resolved experimentally MM- Therefore, 7 h and Jq can be 
measured Il53l . so it would be possible to verify that a finite 
charge current I flows when the condition Ji„ = 0 is met. It 
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is already possible to manipulate locally the temperature of the 
reservoirs IIT8l420l l50ll . thus getting a knob for controlling the 
heat currents and ensuring 7in = 0. 


6. Thermalization induced by a probe reservoir 


Our quantum dot model differs from that of a conventional 
thermocouple in a number of ways; such as the fact that dot 
M has discrete states and that it is in the Coulomb blockade 
regime. Yet, we claim that the exotic power generation is the 
consequence of the non-thermal nature of the distribution in dot 
M, rather than one of the other properties of the dot. 

To verify this, we model the thermalization processes within 
dot M phenomenologically, by considering an additional probe 
reservoir (P), which is tunnel-coupled to dot M, as indicated 
in Fig. The working principle of the probe ll54H5^ is anal¬ 
ogous to the one of a voltage probe ll57ll . with the additional 
feature of probing the temperature. In other words, the probe 
is treated as an electronic reservoir in thermal equilibrium at a 
temperature Tp and a chemical potential yUp, which have to be 
self-consistently determined by requiring that both the charge 
and heat currents, Eq. 0, into the probe reservoir are equal to 
zero, /p = 7p = 0. 

In order to evaluate the currents into the probe and to extract 
Tp and yUp that make them vanish, we treat the probe reservoir 
on the same footing as the other electronic reservoirs introduced 
in Sec.|^ The coupling between dot M and the probe P is 

Ze = TP /p ') , (20) 

with Eq. 0’ s sum over a extended to include P. Then, Ip and 
Jp are given by Eqs. 0 with a = P. 

If the probe is much more weakly coupled to dot M than the 
other two electronic reservoirs, (yp « 7l,r), it does not influ¬ 
ence the dynamics of the 3-dots system, i.e. it does not alter the 
results discussed in Sec.[^and|^ In contrast, if yp is comparable 
with yu.R, or even larger, then the probe mimics phenomenolog¬ 
ically the process of thermalizing transitions within dot M (i.e. 
transitions that take the dot towards thermal equilibrium, with¬ 
out changing its average energy or charge). In other words, the 
state of dot M becomes closer to thermal as we increase the 
coupling to the probe reservoir. 

To measure how far dot M is from a thermal state, we first 
note that if it were in a thermal state corresponding to equilib¬ 
rium with the probe reservoir, then fluctuations of m induced by 
the probe reservoir in both directions (1 —> 0 and 0 —> 1) would 
have the same rate, irrespective of the state of the other dots. 
In other words, dot M is in a thermal state in equilibrium with 
the probe reservoir P if [Tp]]^i^^pihc = [^vtxhcP^hc for all h, c. A 
little algrebra shows that this thermal state corresponds to 


Pi he 

POhe 


exp 


Eihc - Eohe - Pp 




BJ p 


( 21 ) 


for all h,c. Thus, we define the following measure of how far 
dot M is from thermal equilibrium. 


Dfic — 


Plhc 

POhe 


exp 


Elhc ^Ohe Pp 


kpTp 


- 1 


( 22 ) 


Pig. I^a) shows how increasing the coupling yp to the probe 
tends to reduce Dhc for all h, c, indicating that the steady state of 
dot M becomes closer to that of the thermal state with tempera¬ 
ture Tp. Pig.j^b) shows that as dot M’s state becomes closer to 
the thermal state, the power generated at - 0 goes to zero. 
Thus the exotic power generation that we discuss in this arti¬ 
cle is a consequence of the fact that dot M is maintained in a 
non-thermal state by its couplings to reservoirs H and C. 


7. Comparing with a conventional thermocouple 


A conventional thermocouple corresponds to the system in 
Pig. [^a), with the quantum dot replaced by a macroscopic 
reservoir M that is larger than the electron’s thermalization 
length, so it equilibrates at a temperature Tm between Th and 
Tq. The role of the coupling to reservoirs H and C is to establish 
the temperature I'm- In the steady-state. I'm is the temperature 
which ensures that the heat-flow out of reservoir M into the 
thermoelectrics is equal to the heat flow into reservoir M from 
reservoirs H and C, which we call Jm- Given 7in and I'm, the 
thermocouple’s efficiency, rj = Pgen/An, cannot exceed Carnot 
efficiency T/carnot = (1 - TqITm)- Tfius 


^gen — ^Camot ■fir 


(23) 


with the equals sign holding only if the two thermoelectrics are 
Carnot efficient, corresponding to tfiem botfi having a figure of 
merit ZT —> oo. 

We can immediately see that a system like ours, which gen¬ 
erates finite Pgen when Jin = 0, violates the bound in Eq. ( [2^ . 
Here, we wish to see over what range of parameters our system 
violates this bound. To do this, we need to associate a temper¬ 
ature, I'm, to the non-thermal state of the quantum dot, so we 
can calculate T/carnot- There is no unique way of defining suefi 
a temperature, however a natural way to do so is to say it is 
the temperature that a probe would measure if placed in ther¬ 
mal contact with the dot Il54l45^ . As such, we can use probe 
reservoir P considered in the previous section in the limit of 
very weak coupling to the dot (yp « yr.R), so it has negligible 
effect on dot-M’s state. Although dot M is in general in a non- 
thermal state, we can associate to it an effective temperature 
Tm, by defining Tm as the temperature of the probe reservoir 
when there is no charge or heat current flow between dot M and 
reservoir P, i.e. Tm = Tp|/p=yp=o- Note that Ref. Il5^ argued that 
for having a meaningful definition of it is crucial that the 
probe’s coupling to the dot is the same at all energies, so the 
probe is sensitive to the whole energy distribution of the dot. 

Having defined an effective temperature for dot M, and 
knowing the heat-flowing into it, Jin, we can plot the upper 
bound given by Eq. (j2^ for any given parameters of the model. 
This is represented as a dashed curve in all plots, see Eigs. 2]|6 


We see that the power generated by our machine (solid curves) 
exceeds this classical bound over a broad range of parameters, 
including when Tq is equal to or above Tq. 

Note that //camot-fin will be negative whenever Jin and Tm-Tq 
are of opposite sign. Eor example, in Eig.|^one has 7in < 0 for 
TcITo < 0.94, while I'm < Tq only for TqITq < 0.11. Thus 
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(a) How probe couplings affect 

relative oecupations 



E\hc Eofic 


(b) Power generated for increasing (c) Difference between power generated and 

probe coupling upper bound on conventional thermocouple 




Figure 6: The same system parameters as in Fig.[^ but with finite coupling to the temperature probe, (a) How dot-M’s state becomes increasingly thermal as the 
probe coupling is increased. We show a scatter plot of — as defined in Eq. H) — versus (E[hc - -Eo/jc) with {h,c) = (0,0),(1,0),(0,1) and (1,1). As the 
coupling to the probe yp is increased, the scatter of the points reduces towards zero. Here we show yp/y = 0.1,0.5,1,10 (filled diamonds going from black to 
green), with the arrows indicating increasing probe coupling. We also show the two limits; yp 0 (open black diamond, for this we actually take yp = O.OOly) 
and yprobe ^ (filled red diamond, for this we actually take yp = lOOOy). (b) The power output as a function of Tq for different probe couplings (solid curves) 
for yp/y = 0.1,0.5,1,10. The filled circles on each curve indicate the value of Pgen when Jin = 0, we see that this goes towards zero as the coupling to the 
probe increases. The upper bound on conventional thermocouples, t/camot^fin’ indicated by the dashed curves for weak and strong probe coupling (black for 
yp = O.ly and green for yp = lOy). Note that for strong probe coupling the dot power output is below the bound for a conventional thermocouple for all Tq- This is 
seen more clearly in (c), which plots the difference in power generation between our quantum dot and the upper bound on a conventional thermocouple (again for 
yp/y = 0.1,0.5,1,10). 


for 0.77 < TcITo < 0.94, the dot M is hotter than reservoirs L 
and R but, since Jin < 0, we know that energy is flowing from 
reservoirs L and R into the dot. In a macroscopic thermocouple, 
this can occur only if one provides energy to the system by 
replacing the load with a power supply that forces the current 
flowing through the circuit. In this case, power is not generated, 
but absorbed by the system, i.e. Pgen < 0. In contrast, our 
quantum dot system generates power, Pgen > 0, even in this 
range of parameters, for the reasons discussed in Sec. 

Figs.|^b) and (c) show that as dot M’s state becomes closer 
to that thermal state, the violation of Eq. (231 goes away, leaving 
a system with Pgen = rjJm, where rj is the system’s efficiency, 
and it is less than the Carnot efficiency ?7camot- For example, for 
7 p = lOy in Fig. |^b), tj is of order O.25?7camot- This confirms 
our conclusion that a non-thermal steady-state distribution of 
electrons in dot M is necessary for it to generate more power 
than any conventional thermoelectric could. 


8. Conclusions 

We have investigated the unexpected properties that arise in 
thermoelectric conductors when they become smaller than the 
lengthscale in which electrons thermalise. For this, we have 
considered a quantum dot based thermocouple in a four termi¬ 
nal geometry; two terminals support the electrical current, with 
the other two being the hot and cold environments with which 
energy is exchanged. Our system can supply electrical power 
to a load. Surprisingly, this can be done also in configurations 
in which no net heat is absorbed from the thermal reservoirs. 
This exotic power generation can be explained in terms of non¬ 
local Coulomb heat drag in which the heat flow from the hot 
to the cold reservoir induces both a heat and a charge flow in 
the capacitively coupled electronic circuit. This effect should 
be realizable in systems very close to those in recent experi¬ 
ments GlHlol, as outlined in Sec. |5.3| 


We argue that this exotic power generation relies on the fact 
that quantum dot M is maintained in an non-thermal state, by 
showing that it disappears when relaxation effects destroy the 
non-thermal state. As such it cannot occur in a conventional 
macroscopic thermocouple heat-engine. We compare our quan¬ 
tum dot thermocouple to its macroscopic equivalent, and And 
that avoiding the thermalization of carriers can help to improve 
the thermoelectric performance. In particular, the generated 
power can be larger than that of any macroscopic thermocou¬ 
ple, even when the latter works with Carnot efficiency. In the 
long term, we hope that this observation will be applied in ther¬ 
moelectric and photovoltaic applications, since in both cases 
a large part of the heat flowing into the central region of the 
thermocouple (or the central region of the p-n junction in pho- 
tovoltaics) flows out again into the cold environment, without 
contributing to the power generation. It would be a great ben¬ 
efit if a part of this “lost” energy flow could be used to make 
electrical power in the manner presented here. 

The same effect could be achieved with non-interacting con¬ 
ductors in magnetic fields, via a Nernst effect lIS^IMI . Note that 
Refs. 15911^ assumed that the reservoirs that carry the charge 
current (reservoirs L and R in our Fig. [T]| are at temperatures 
chosen such that there is no heat flow out of them. In this situ¬ 
ation, energy conservation guarantees that Pgen = 0 for Ji„ = 0. 
However, for other temperatures, one expects to observe the ex¬ 
otic power generation that is discussed in the present paper. 

An open question is how to engineer the exotic power gen¬ 
eration when the capacitative coupling is replaced by exchange 
of photons ll58ll or phonons Il9l [3Tl - l3^ . Not only is this crucial 
for thermoelectric and photovoltaic applications, it would help 
to elucidate the fundamental requirements for the effect. 

As a step in this direction, it may be worthwhile studying 
the linear-response regime in more detail. The results pre¬ 
sented here are valid for arbitrary Ph, Tq and Tq-, while the 
observed effect is most prominent for very different tempera- 
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tures, it also exists when the temperature differences are small. 
In this limit, one could analyze the effect using a multi-terminal 
Onsager treatment MM, which is more universal, but only 
works when there is a linear relationship between thermody¬ 
namic forces and currents. 
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